关于求解随机用户均衡问题的截断拟牛顿型信赖域法研究

信赖域法是一种保证全局收敛性的优化算法,为避免Hessian矩阵的计算,基于拟牛顿校正公式构造了求解带线性等式约束的非线性规划问题的截断拟牛顿型信赖域法.首先给出了截断拟牛顿型信赖域法的构造过程及具体步骤;然后针对随机用户均衡模型中变量和约束的特点对算法进行了修正,并将多种拟牛顿校正公式下所得结果与牛顿型信赖域法的结果进行了比较,结果发现基于对称秩1校正公式的信赖域法更为合适.最后基于数值算例结果得到了一些在算法编程过程中的重要结论,对其它形式信赖域法的编程实现具有一定的参考意义.     

红外无损探测中多宗量多热源反演问题的研究

为研究红外无损探测稳态多热源反演逆问题, 建立不同形状的均质与非均质稳态热传导模型, 其中内热源个数、位置、强度、面积均为未知项. 基于数值算法中有限元算法对模型进行离散分析, 化简有限元矩阵方程, 最终转化为对Ax = b高度欠定方程的求解. 首次利用分段多项式谱截断奇异值分解法处理内热源逆问题, 并对算法进行改进, 有效改善了该算法在处理多热源反演时存在的严重的热源叠加效应. 根据反演出的内热源信息, 利用有限元算法计算重构出整个模型内所有节点的温度分布. 运用数值仿真Comsol软件和具体实物实验对算法进行有效性评估, 并验证算法在不同热传导模型中的表现. 结果表明, 算法能够准确反演出多热源各参量信息, 在非均质材料模型中仍能准确地反演出热源项, 并有效重构出模型内温度场. 该算法可应用于材料无损检测及人体红外医学成像等领域.     

Stress-strain analysis of a composite truncated conical shell

The stress-strain state of a nonshallow truncated conical shell made of a composite material and subjected to axial forces at the ends is studied. A variational difference method is used. Calculations are carried out for an orthotropic shell with low shear stiffness __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 7, pp. 66–72, July 2007.     

Chaotic Oscillation of Saddle Form Cable-Suspended Roofs under Vertical Excitation Action

The purpose of this paper is to investigate the dynamic behaviour of saddle form cable-suspended roofs under vertical excitation action. The governing equations of this problem are system of nonlinear partial differential and integral equations. We first establish a spectral equation, and then consider a model with one coefficient, i.e., a perturbed Duffing equation. The analytical solution is derived for the Duffing equation. Successive approximation solutions can be obtained in likely way for each time to only one new unknown function of time. Numerical results are given for our analytical solution. By using the Melnikov method, it is shown that the spectral system has chaotic solutions and subharmonic solutions under determined parametric conditions.     

NEW METHOD FOR LOW ORDER SPECTRAL MODEL AND ITS APPLICATION

In order to overcome the deficiency in classical method of low order spectral model, a new method for low order spectral model was advanced. Through calculating the multiple correlation coefficients between combinations of different functions and the recorded data under the least square criterion, the truncated functions which can mostly reflect the studied physical phenomenon were objectively distilled from these data. The new method overcomes the deficiency of artificially selecting the truncated functions in the classical low order spectral model. The new method being applied to study the inter-annual variation of summer atmospheric circulation over Northern Hemisphere, the truncated functions were obtained with the atmospheric circulation data of June 1994 and June 1998. The mechanisms for the two-summer atmospheric circulation variations over Northern Hemisphere were obtained with two-layer quasi-geostrophic baroclinic equation.     

THE CONSTRUCTION OF WAVELET-BASED TRUNCATED CONICAL SHELL ELEMENT USING B-SPLINE WAVELET ON THE INTERVAL

Based on B-spline wavelet on the interval (BSWI), two classes of truncated conicalshell elements were constructed to solve axisymmetric problems, i.e. BSWI thin truncated conicalshell element and BSWI moderately thick truncated conical shell element with independent slope-deformation interpolation. In the construction of wavelet-based element, instead of traditionalpolynomial interpolation, the scaling functions of BSWI were employed to form the shape functionsthrough the constructed elemental transformation matrix,and then construct BSWI element viathe variational principle. Unlike the process of direct wavelets adding in the wavelet Galerkinmethod, the elemental displacement field represented by the coefficients of wavelets expansionwas transformed into edges and internal modes via the constructed transformation matrix. BSWIelement combines the accuracy of B-spline function approximation and various wavelet-basedelements for structural analysis. Some static and dynamic numerical examples of conical shellswere studied to demonstrate the present element with higher efficiency and precision than thetraditional element.     

Moving boundary truncated grid method: Multidimensional quantum dynamics

The moving boundary truncated grid (TG) method is used to study wave packet dynamics of multidimensional quantum systems. As time evolves, appropriate Eulerian grid points required for propagating a wave packet are activated and deactivated with no advance information about the dynamics. This method is applied to the Henon-Heiles potential and wave packet barrier scattering in two, three, and four dimensions. Computational results demonstrate that the TG method not only leads to a great reduction in the number of grid points needed to perform accurate calculations but also is computationally more efficient than the full grid calculations.     

基于COMSOL软件下截锥体电阻的计算

本文应用Comsol软件讨论了截锥体内电势和电流的分布,并计算了截锥体的电阻.     

删失截断情形下失效率变点模型的Bayes参数估计

通过添加部分缺失寿命变量数据,得到了删失截断情形下失效率变点模型相对简单的似然函数.讨论了所添加缺失数据变量的概率分布和随机抽样方法.利用Monte Carlo EM算法对未知参数进行了迭代.结合Metropolis-Hastings算法对参数的满条件分布进行了Gibbs抽样,基于Gibbs样本对参数进行估计,详细介绍了MCMC方法的实施步骤.随机模拟试验的结果表明各参数Bayes估计的精度较高.     

A Mixture Model of Truncated Zeta Distributions with Applications to Scientific Collaboration Networks

The degree distribution has attracted considerable attention from network scientists in the last few decades to have knowledge of the topological structure of networks. It is widely acknowledged that many real networks have power-law degree distributions. However, the deviation from such a behavior often appears when the range of degrees is small. Even worse, the conventional employment of the continuous power-law distribution usually causes an inaccurate inference as the degree should be discrete-valued. To remedy these obstacles, we propose a finite mixture model of truncated zeta distributions for a broad range of degrees that disobeys a power-law behavior in the range of small degrees while maintaining the scale-free behavior. The maximum likelihood algorithm alongside the model selection method is presented to estimate model parameters and the number of mixture components. The validity of the suggested algorithm is evidenced by Monte Carlo simulations. We apply our method to five disciplines of scientific collaboration networks with remarkable interpretations. The proposed model outperforms the other alternatives in terms of the goodness-of-fit.